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Abstract 

We calculate the next-to-leading QCD corrections to the effective Hamiltonian for B — > 
Xse~^e~ in the NDR and HV schemes. We give for the first time analytic expressions for 
the Wilson Coefficient of the operator Qg = {sb)v-A{ee)v in the NDR and HV schemes. 
Calculating the relevant matrix elements of local operators in the spectator model we 
demonstrate the scheme independence of the resulting short distance contribution to 
the physical amplitude. Keeping consistently only leading and next-to-leading terms, 
we find an analytic formula for the differential dilepton invariant mass distribution in 
the spectator model. Numerical analysis of the mt, and 0{m\^) dependences 
of this formula is presented. We compare our results with those given in the literature. 



* Supported by the German Bundesministerium fiir Forschung und Technologie under contract 06 TM 
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1 Introduction 



The rare decay B Xse~^e~ has been the subject of many theoretical studies in the frame- 
work of the standard model and its extensions such as the two Higgs doublet models and 
models involving supersymmetry In particular the strong dependence of B ^ X^e^e" 

on rrit has been stressed by Hou et. al. It is clear that once B — >■ X^e^e" has been ob- 
served, it will offer an useful test of the standard model and of its extensions. To this end the 
relevant branching ratio, the dilepton invariant mass distribution and other distributions of 
interest should be calculated with sufficient precision. In particular the QCD effects should 
be properly taken into account. 

The central element in any analysis of -B — * Xse~^e~ is the effective Hamiltonian for 
AB = 1 decays relevant for scales fi ~ (9(mb) in which the short distance QCD effects 
are taken into account in the framework of a renormalization group improved perturbation 
theory. These short distance QCD effects have been calculated over the last years with 
increasing precision by several groups [@, P,[lO| culminating in a complete next-to- leading 
QCD calculation presented by Misiak in ref. W^ and very recently in a corrected version 



in H. 

The actual calculation of B ^ X^e+e" involves not only the evaluation of Wilson coeffi- 
cients of ten local operators (see which mix under renormalization but also the cal- 
culation of the corresponding matrix elements of these operators relevant for B Xse^e~ . 
The latter part of the analysis can be done in the spectator model, which, according to 
heavy quark effective theory, for B-decays should offer a good approximation to QCD. One 
can also include the non-perturbative (9(l/m^) corrections to the spectator model which en- 



hance the rate for B X^e^e" by roughly 10% [|l^]. A realistic phenomenological analysis 
should also include the long distance contributions which are mainly due to the J/ ip and ip' 
resonances |]14|-[T6[]. Since in this paper we are mainly interested in the next-to- leading short 
distance QCD corrections to the spectator model we will not include these complications in 
what follows. 

It is well known that the Wilson coefficients of local operators depend beyond the leading 
logarithmic approximation on the renormalization scheme for operators, in particular on 
the treatment of 75 in D 7^ 4 dimensions. This dependence must be cancelled by the 
scheme dependence present in the matrix elements of operators so that the final decay 
amplitude does not depend on the renormalization scheme. In the context of B ^ Xse^e~ 
this point has been emphasized in particular by Grinstein et. al. 0. Other examples such 
as, K ^ 7T7C, Kl s ~* vr°e^e^, B Xs'j can be found in refs. [17-19|. The interesting 
feature of B XgC^e^ as compared to decays such as K vrvr, is the fact that due to the 
ability of calculating reliably the matrix elements of all operators contributing to this decay, 
the cancellation of scheme dependence can be demonstrated in the actual calculation of the 
short distance part of the physical amplitude. 

Now all the existing calculations of i? ^ X^e+e" use the NDR renormalization scheme 
(anticommuting 75 in D 7^ 4 dimensions). Even if arguments have been given, in particular 
in and how the cancellation of the scheme dependence in B —>■ Xse~^e~ would 
take place, it is of interest to see this exphcitly by calculating this decay in two different 
renormalization schemes. In addition, in view of the complexity of next-to-leading order 
(NLO) calculations and the fact that the only complete NLO analysis of -B ^ X^e^e^ has 
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been done by a single person, it is important to check the results of refs. [|n|,|T^. 

Here we will present the calculations of the Wilson coefficients and matrix elements relevant 
for B — > XsC+e" in two renormalization schemes (NDR and HV demonstrating the 
scheme independence of the resulting amplitude. Beside this the main results of our paper 
are as follows: 

• We give for the first time analytic NLO expressions for the Wilson coefficient of the 
operator Qg = (s6)y_yi(ee)v in the NDR and HV schemes. 

• Calculating the matrix elements of local operators in the spectator model we fully 
agree with Misiak's result for the dilepton invariant mass distribution very recently 
given in [O. 



We find, that in the HV scheme the scheme dependent term in the matrix elements (the 
so called ^-term) receives in addition to current-current contributions also contributions 
from QCD penguin operators which are necessary for the cancellation of the scheme 
dependence in the final amplitude. This should be compared with the discussion 



of the scheme dependence given in refs. and [11| where the ^-term received only 
contributions from current- current operators. 

We stress that in a consistent NLO analysis of the decay B X^e^e" , one should 
on one hand calculate the Wilson coefficient of the operator Qg = (s6)v-A(ee)v' in- 
cluding leading and next-to-leading logarithms, but on the other hand only leading 
logarithms should be kept in the remaining Wilson coefficients. Only then a scheme 
independent amplitude can be obtained. This special treatment of Qg is related to the 
fact that strictly speaking in the leading logarithmic approximation only this operator 
contributes to i? — > X^e+e" . The contributions of the usual current-current opera- 
tors, QCD penguin operators, magnetic penguin operators and of Qio = (s6)v-A(ee)A 
enter only at the NLO level and to be consistent only the leading contributions to the 
corresponding Wilson coefficients should be included. In this respect we differ from the 



original analysis of Misiak [|TT| who in his numerical evaluation of B ^ Xse^e" also 
included partially known NLO corrections to Wilson coefficients of operators Qi[i 7^ 9). 
These additional corrections are, however, scheme dependent and are really a part of 
still higher order in the renormalization group improved perturbation theory. The most 



recent analysis of Misiak [|T2[ does not include these contributions and can be directly 



compared with the present paper. 

• Keeping consistently only the leading and next-to-leading contributions to 5 ^ 
Xse^e~ we are able to give analytic expressions for all Wilson coefficients which should 
be useful for phenomenological applications. 

Our paper is organized as follows: 

In sect. 1^ we collect the master formulae for B X^e+e^ in the spectator model which 
include consistently leading and next-to- leading logrithms. In sect. | we describe some 
details of the NLO calculation of the Wilson coefficient Cg{fi) and of the relevant one- loop 
matrix elements in NDR and HV schemes. In sect. ^ we present a numerical analysis. We 
end our paper with a brief summary of the main results. 
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2 Master Formulae 



2.1 Operators 



(2.1) 



Our basis of operators is given as follows: 

Ql = {SaCf3)v~A{C/3bQ)v-A 

Q2 = isc)v-Aicb)v-A 

Qs = isb)v-AJ2qiqq)v~A 

Qi = {Sabp)v~AY.q{<if3<la)v-A 

Q5 = {sh)v~-AY.q{qq)v+A 

Qe = iSabi3)v-Aj2q{qi3qa)v+A 

Qg = {sb)v-A{ee)v 
Qw = {sb)v-A{ee)A 

where a and f3 denote colour indices. We omit the colour indices for the colour-singlet 
currents. Labels (V ± A) refer to 7^(1 ± 75). (5i,2 are the current- current operators, Qs-e 
the QCD penguin operators, Qi^s "magnetic penguin" operators and (^9,10 semi-leptonic 
electroweak penguin operators. Our normalizations are as in refs. and . 



2.2 Wilson Coefficients 

The Wilson coefficients for the operators Q1-Q7 are given in the leading logarithmic approx- 



imation by |[T8|, |2lH23[ 

Cfifi) = jzk.^r,^^ (J = l,...6) 



with 



1=1 



i=l 



V 

Cf\M^) 



(2.2) 
(2.3) 



(2.4) 

(2.5) 
(2.6) 
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where Xt = m'^/M^ and A{x) and F{x) are defined in ( 2.14 ) and ( 2.19 ). The numbers a, 
kji and hi are given by 



tti = ( 


14 


16 


_6_ 

9"? ' 


_ 12 
9*? ' 


0.4086, 


-0.4230, 


-0.8994, 


0.1456 


hi = ( 


0, 


0, 


1 

2 ' 


1 

2 ' 


0, 


0, 


0, 





A;2i = ( 


0, 


0, 


1 

2' 


1 

2' 


0, 


0, 


0, 





hi = ( 


0, 


0, 


1 

14' 


1 

6' 


0.0510, 


-0.1403, 


-0.0113, 


0.0054 


= ( 


0, 


0, 


1 

14' 


1 

6' 


0.0984, 


0.1214, 


0.0156, 


0.0026 


= ( 


0, 


0, 


0, 


0, 


-0.0397, 


0.0117, 


-0.0025, 


0.0304 


hi = ( 


0, 


0, 


0, 


0, 


0.0335, 


0.0239, 


-0.0462, 


-0.0112 


/ij = ( 


2.2996, 


-1.0880, 


3 
7' 


1 

14' 


-0.6494, 


-0.0380, 


-0.0186, 


-0.0057 



(2.7) 

The first correct calculation of the two-loop anomalous dimensions relevant for ( |2.3| ) has 



been presented in pl|,p^ and confirmed subsequently in |[T^, p4| , |25[| . 



The coefficient Cg'^'^^^{fi) does not enter the formula for B —>■ X^e'^e at this level of 



accuracy. An analytic formula is given in ref. 
The coefficient of Qio is given by 



Cio(M, 



a 



w 



7rCio(Mw), 



Cio(M, 



Y(xt) 



wj 



sin 9^^ 



(2. 



with Y{x) given in (|2.13| ). Since Qio does not renormahze under QCD, its coefficient does not 
depend on ^ ^ 0{m^). The only renormalization scale dependence in ( |2.8| ) enters through 
the definition of the top quark mass. We will return to this issue in sect. |. 
Finally, including leading as well as next-to-leading logarithms, we find 



with 



r^NDR 
^^9 



r<NDR 



a 
2^ 



Po^^^ + - 4Z(x,) + P^E{x,) 

sm Biy 



(2.9) 
(2.10) 



Here 



NDR 




Pe 
Y(x) 



-0.1875 + ^Pir7"'+^; 



i=l 



asiMw) 
+ 1.2468 + V"' [r^^^ + SiV] 



i=l 



0.1405 + ^gir7"'+^ 



i=l 



C{x)-B{x), Z{x)=C{x) + -D{x). 



A{x) 



x{8x^ + 5x-7) x^{2 - 3x) 
12(x-l)3 + 2(x-l)4 ' 



(2.11) 

(2.12) 
(2.13) 

(2.14) 
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C{x 
D{x 
E{x 
Fix 



X 



+ 



X 



4(1 -x) 4(a;-l)2 
x{x — 6) x{3x + 2) 



\nx, 



(x-l) 8(a;-l)2 
-19x3 + 25x2 x2(5x2-2x-6) 



H ^ — ; r- — - In X In X, 

18 x-l 4 9 



36(x- 1)3 
x(18-llx-x2) , x2(15 - 16x + 4x2 

12(1 -x)3 
xfx^ — 5x — 21 



+ 



6 1 - X 



In X In X, 

3 



4(x- 1)3 



+ 



3x^ 



2(x - 11 



Inx. 



(2.15) 
(2.16) 
(2.17) 
(2.18) 
(2.19) 



The coefficients Pi, rf^^, Sj, and are found to be as follows: 



Pi = ( 


0, 


0, 


80 

203' 


8 

33' 


0.0433, 


0.1384, 


0.1648 


-0.0073 




0, 


0, 


0.8966, 


-0.1960, 


-0.2011, 


0.1328, 


-0.0292, 


-0.1858 


Si = ( 


0, 


0, 


-0.2009, 


-0.3579, 


0.0490, 


-0.3616, 


-0.3554, 


0.0072 


= ( 


0, 


0, 


0, 


0, 


0.0318, 


0.0918, 


-0.2700, 


0.0059 



(2.20) 

Pe is O(10~2) and consequently the last term in (|2.10| ) can be neglected. We keep it 
however in our numerical analysis. 
In the HV scheme only the coefficients are changed. They are given by 

rf^ = ( 0, 0, -0.1193, 0.1003, -0.0473, 0.2323, -0.0133, -0.1799 ). 

Equivalently we can write 



with 



We note that 



pHv ^ pNDR ^ ^Hv:t (3c(o) + - - sa 



(0) 



0, NDR 
-1, HV. 



= 0.1875 

i=l 
8 



= -0.1405, 



i=l 



-1.2468 + -(1 + 0, 



i=l 



i=l 



11 = - 



16 
69' 



(2.21) 

(2.22) 
(2.23) 

(2.24) 
(2.25) 



In this way for rj = 1 we find Pe = 0, P^^^ = 4/9 and P^^ = in accordance with the 
initial conditions in (|3.3|) . Moreover, the second relation in (|2.25| ) assures the correct large 
logarithm in Pq^^^, i- 8/9 ln(Mw//u). The derivation of (|2]^)-(|2]^) is given in sect. |. 



2.3 The Differential Decay Rate 

Introducing 

. _ {Pe+ +Pe-f 
^ 2 

mt 



rric 



(2.26) 
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and calculating the one-loop matrix elements of Qi using the spectator model in the NDR 
scheme we find 



se^e 



where 



Here 



with 



ceu 



a 

47]-2 



Vt 



ts 



2 



4 1 + 



f{z)K{z) 
|(^(0)e//|2 



;i + 2s)(|G 



e//|2 
9 I 



10 1 



12Cf^^^ Re a" 



CI" 



C!:^DRfj{s) + h{z, s) (3Ci°) + ) + 3CP + Cf^ + 3Cf ) 



^9 

I l,\ 1 . I 1 "H-l . 



(0) 



(0) 



3d°) + 



8^ rub 8, 8 4 

-- In H 1 — X 

9 u 9 27 9 



[2 + x)|l 



1 1/2 



(in 



ITC 



2 arctan ■ 



-1' 



for X 
for x 



s 
4^2 



< 1 
>1, 



h{0,s) 
?7(s) 



8 8, mb 4, . 4. 
In In s H — tiT. 



27 9 /i 

l-8z^ + 8z^ 
2a^(/i) 



1 - 

1 + 



9 

- 24z^ In z, 



TT 



--TT 

9 



- In s ln(l 



2s(l + s)(l 



3(1 - s)2(l + 2s) 



2s) 5 
m s + 



-s)- 
9s 



5 + 4s 
3(1 + 2s) 



ln(l - s) 



6(l-s)(l + 2s) 



+ 



(2.27) 



(2.28) 



(2.29) 



(2.30) 
(2.31) 
(2.32) 

(2.33) 



(2.34) 



Here f{z) and k,{z) are the phase-space factor and the single gluon QCD correction to the 
b ceu decay |^,^ respectively, fj on the other hand represents single gluon corrections 
to the matrix element of Qg with mg = ]12|,^. For consistency reasons this correction 
should only multiply the leading logarithmic term in P^^^. 

In the HV scheme the one-loop matrix elements are different and one finds an additional 
explicit contribution to ( |2.28D given by 



HV 



9 



3Cr + a' 



(0) 



3d°) 



(2.35) 



However C^^^ has to be replaced by C|^^ given in ( |2.1(J| ) and ( |2.22| ) and consequently Cg'^^ 
is the same in both schemes. 
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The first term in the function h{z, s) in ( p^.29| ) represents the leading /z-dependence in the 
matrix elements. It is cancelled by the /x-dependence present in the leading logarithm in Cg. 
The /x-dependence present in the coefficients of the other operators can only be cancelled by 
going to still higher order in the renormalization group improved perturbation theory. To 
this end the matrix elements of four-quark operators should be evaluated at two-loop level. 
Also certain unknown three-loop anomalous dimensions should be included in the evaluation 
of Cj^^ and Cg [|18],[19|. Certainly this is beyond the scope of the present paper and we will 



only investigate the left-over //-dependence in sect. |[ 

The fact that the coefficient Cg should include next-to-leading logarithms and the other 
coefficients should be calculated in the leading logarithmic approximation is easy to un- 
derstand. There is a large logarithm in Cg represented by l/a^ in Pq in Conse- 
quently the renormalization group improved perturbation theory for Cg has the structure 
0{l/as) + 0{1) + 0{as) + . . . whereas the corresponding series for the remaining coeffi- 
cients is 0{1) + 0{as) + . . .. Therefore in order to find the next-to-leading 0{1) term, the 
full two-loop renormalization group analysis for the operators in ( |2.1|) has to be performed 
in order to find Cg, but the coefficients of the remaining operators should be taken in the 
leading logarithmic approximation. This is gratifying because the coefficient of the mag- 
netic operator Qi is known only in the leading logarithmic approximation. Qi does not mix 
with Qg and has no impact on the coefficients Ci-Cq. Consequently the necessary two- loop 
renormalization group analysis of Cg can be performed independently of the presence of the 
magnetic operators, which was also the case of the decay — > 7r°e"'"e~ presented in ref. 

Let us finally compare our main formulae (|2.27|) -( ^.35|) with the ones given in the literature: 



i) The general expression ( 2.27|) with k,{z) = 1 is due to Grinstein et. al. who in 



their approximate leading order renormalization group analysis kept only the operators 
Qi, Q2, Q7, Q9, Qio- 



ii) Inserting cf ^ and C^^^ in Qand into { ^71^ we find an analytic expression 



for Cg^^ which agrees with a recent independent calculation of Misiak |jT2[ 



iii) The sign of iir in ( p^.29D differs from the one given in ^ and [|rT|] but agrees with []T^ 



and also with the work of Fleischer 29 . 



iv) The "^-term" given in ( p.35| ) contains in the HV scheme also contributions from the 
operators Q3 and Q4, which are however negligible. The discussion of the "^-term" in 
refs. and [|ri|] does not apply then to the HV scheme. 



3 Technical Details 



3.1 Wilson Coefficients 



In order to calculate the coefficient Cg including next-to-leading order corrections we have 
to perform in principle a two-loop renormalization group analysis for the full set of operators 



given in However, Qio is not renormalized and the dimension five operators Qj and 

Qs have no impact on Cg. Consequently only a set of seven operators, Qi-q and Qg, has to 
be considered. This is precisely the case of the decay — > n^e^e~ considered in [^] except 
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for an appropriate change of quark flavours and the fact that 0(mb) instead of 

fj, ~ C(lGeV) should be considered. Because our detailed NLO analysis of Kl 7r°e"'"e~ 
has already been published we will only discuss very briefly an analogous calculation of 
B Xse~^e~ , referring the interested reader to ||19||. We should stress that Misiak [jll],|12 
used different conventions for the evanescent operators than used in and here. The 
agreement on Cg^^ is therefore particularly satisfying. 

Integrating out simultaneously W, Z and t we construct first the effective Hamiltonian for 
AB = 1 transitions relevant for b se^e~ with the operators normalized at = Mw. 
Dropping the operators Q?, Q% and Qio for the reasons stated above and using the unitarity 
of the CKM matrix we find 



K//(A5 = 1) 



+ 



Gp 
G 



J2G,{M^)Q, + G'g{MM 



.1=1 



[Ci(Mw)(Q1"^ - Qi) + C2(Mw)(Q5"^ - Q2)] ■ (3.1 



Here Q^^l are obtained from Q12 through the replacement c u. In order to make all the el- 
ements of the anomalous dimension matrix be of the same order in a^, we have appropriately 
rescaled Cg and Qg: 

Q'g = ^Qg, CM = ^Cgi^i) (3.2) 

Note that because of GIM cancellation there are no penguin contributions in the term propor- 
tional to V*gVub- They would appear only at scales fi < rric as was the case in —>■ n^e'^e~ . 
Since |14*sV^6/V^sVtb| < 0.02 we will drop the second term in what follows. 

The initial conditions at /i = Mw for the coefficients Ci-Cq in NDR and HV schemes have 
been given in sect. 2.4 and in the appendix A of ref. []T^ respectively. Here it suffices to give 
only the initial condition for the coefficient C'g (denoted by C'^y in |T^) which reads: 



C^(Mw) 



WJ 



271 



sm y 



(3.3) 



where ^ has been defined in ( p.23|) . The Xt dependence originates in box diagrams and in 
the 7- and Z-penguin diagrams ||30| . 
With 

C^^iG„...,Ge,G'g) (3.4) 

one can calculate the coefficients Cj(/i) by using the evolution operator U^^fijMyj) relevant 
for an effective theory with / = 5 flavours: 



C(/i) = ?75(/i,Mw)C(Mw) 



(3.5) 



An explicit expression for Ur, is given in sect. 2 of [|T9|] where also the relevant expressions for 
one- and two-loop anomalous dimensions can be found. One only has to set / = 5, u = 2 
and d = 3 in the formulae given in [|l^ . 
Using (|3.5D and rescaling back the operator Qg we find at yU ~ C(mb) 



HeffiAB = 1) = -^V*V,, 



^a(/i)g, + C9(/i)Qg 



.4 = 1 



(3.6) 
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with the coefficient Cq{^) given in ( |2.10| ) and ( ^.22| ) for NDR and HV schemes respectively. 
The result for HV can either be found directly using (|3.5|) or by using the relation 



C^y(^^) =ii- ^Af^ J C7^^«(/i) (3.7) 



with the matrix Af given in appendix A of ref. |T9 |. 



3.2 One-Loop Matrix Elements 

The operators Qj and Qio contribute at this level of accuracy only through tree level matrix 
elements. Qs contributes only through the renormalization of and its impact is only felt 
in cf^^^^ . The four-quark operators Qi-e, contribute at one-loop level through the diagrams 
in fig. m where "® ®" denotes the operator insertion. Finally at next-to-leading level 0{as) 
corrections to the matrix element (Qg) have to be calculated. 



a) b) 

Figure 1: The two possibilites for insertion of a four-quark operator into a penguin diagram. 

Let us begin with (Qi-e)- As usual two types of insertions of the operators into the 
penguin diagrams have to be considered. As already discussed in ref. the appearance 
of a closed fermion loop in fig. ||a does not pose any problems in the NDR scheme because 
nowhere in the calculation one has to evaluate Tr[7^7y7p7A75]. The diagrams in fig. |T]have 
been evaluated for the operators Qi and Q2 by Grinstein et. al. and by Misiak [jlT| for the 
full set Qi-Qq. These calculations have been done in the NDR scheme. Calculating these 
diagrams in the NDR and HV schemes we find 

(Qi) = 1^ (sM^, s) - (Q9)o 



(Q2) = 1^ (M^, S) - (Q9)0 



o / 1 2 4 \ 

m = —(3h{z,s)-2h{l,s)--h{0,s) + - + -^j{Q,)o (3. 

(Q4) = ^(h{z,s)-2h{l,s)-h{0,s) + ^ + ^^yQ,)o 

m = ^{3h{z,s)-hil,s) + fj{Q,)o 



(Qe) = ^{hiz,s)-^h{l,s) + '^yQ,), 
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with ^ defined in ( ^.23| ), {Q9)q denoting the tree level matrix element of Qg and 



h{z,s) = ^G{z,s)-^-l\n'^. (3.9) 

Here 

1 

G{z, s) = -4 Jdxx{l - x) In (z^ - sx{l - x)) (3.10) 



with z and s defined in (|2.26|) . 
A few remarks should be made: 

• h{z,s), h{l,s) and h{0,s) correspond to internal c, b and massless {u,d,s) quarks in 
fig. |l] respectively. 

• The contributions of {u, d, s) to diagram [l|a) cancel each other and consequently h{0, s) 
represents the contribution of the internal strange quark in diagram |I|b). 

• We note that (Q^) and (Qq) matrix elements do not contain the ^-term. We should 
however stress that generally it is certainly possible to find schemes in which (Q5) and 
(Qq) matrix elements can differ from the ones given in ( |3.8| ). Similarly we have no 
argument that in schemes different from NDR and HV the matrix elements are found 
simply by changing the value of ^ in the formulae given above. It could be that the 
changes are more involved. Consequently the discussions of the ^-term presented in [0 



and 11] are not generally valid. 



The one gluon correction to the matrix element of Qg, ^7(5), can be inferred from p8[ as 



has been noticed by Misiak in []T2[. In |^ a left-handed current has been considered. Thus 
we rewrite the vector current as a sum of left- and right-handed currents. Neglecting the 
electron masses these two contributions do not interfere. Charge conjugation transforms the 
right-handed current into a left-handed one. Since s is invariant under this transformation 
both currents lead to the same invariant mass spectrum. Therefore we can write 

-2 } ~ 

uj{s) = — dxFi{x, s) (3.11) 



(1_ 5)2(1 + 2s) 

with Fi{x, s) defined explicitly in eq. (3.9) of [Q. Calculating the integral we arrive at the 
result given in (|2.34|) which furthermore agrees with Misiak [5^. 



4 Numerical Analysis 



In our numerical analysis we will use 



a 



4:71 



1 - 



(3, lnln(/iVA|g, 



/5oln(/xVA|3) [ PI ln(/xVAy J 



(4.12) 
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with Pq = 23/3 and Pi = 116/3 as appropriate for five flavours. We also take Aj^ = 
(225 ±85) MeV corresponding to as{Mz) = 0.117 ±0.007. For the remaining parameters we 
take 

a = 1/129, rric = 1.4 GeV, 

sin^Ow = 0.23, = 4.8 GeV, (4.13) 

\Vts/V,b\ = 1, Mw = 80.0 GeV. 

In table |l] we show the constant Pq in ( p.ll|) for different fi and Aj^, in the leading order 
corresponding to the first term in ( ^.11| ) and for the NDR and HV schemes as given by ( |2.11| ) 
and ( |2.22| ) respectively. In table ^we show the corresponding values for C*9(/i). To this end 
we set rrit = 170 GeV. 







= 0.140 GeV 




= 0.225 GeV 




= 0.310 GeV 


fi[GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


2.5 


2.052 


2.927 


2.796 


1.932 


2.845 


2.758 


1.834 


2.774 


2.726 


5.0 


1.851 


2.623 


2.402 


1.787 


2.589 


2.394 


1.735 


2.560 


2.387 


7.5 


1.673 


2.389 


2.125 


1.630 


2.371 


2.126 


1.596 


2.356 


2.126 


10.0 


1.524 


2.202 


1.910 


1.493 


2.192 


1.915 


1.468 


2.183 


1.919 



Table 1: The coefficient Pq of Cg for various values of Aj^ and /i. 







= 0.140 GeV 




= 0.225 GeV 




= 0.310 GeV 


/i[GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


2.5 


2.052 


4.495 


4.364 


1.932 


4.413 


4.326 


1.834 


4.341 


4.293 


5.0 


1.851 


4.193 


3.972 


1.787 


4.159 


3.963 


1.735 


4.130 


3.956 


7.5 


1.673 


3.960 


3.696 


1.630 


3.942 


3.696 


1.596 


3.926 


3.697 


10.0 


1.524 


3.774 


3.482 


1.493 


3.763 


3.486 


1.468 


3.754 


3.490 



Table 2: Wilson coefficient Cg for mt = 170 GeV and various values of Aj^ and /x. 
We observe: 

• The NLO corrections to Pq enhance this constant relatively to the LO result by roughly 
45% and 35% in the NDR and HV schemes respectively. This enhancement is analogous 
to the one found in the case of 7r°e+e~ . 

• In calculating Pq in the LO we have used as{^) at one- loop level. Had we used the two- 
loop expression for ^^(/i) we would find for = 5 GeV and A^ = 225 MeV the value 
Pq'^ f» 1.98. Consequently the NLO corrections would have smaller impact. Ref. 
including the next-to-leading term 4/9 would find Pq values roughly 20% smaller than 
Pq^^ given in tab. [l|. 

• It is tempting to compare Pq in table [1| with that found in the absence of QCD 
corrections. In the limit — > we find Pq^^ = 8/9 ln(Mw//u) ±4/9 and 
pHV = 3/9 in(Mw//i) which for /z = 5 GeV give P^^^ = 2.91 and P^^ = 2.46. Com- 
paring these values with table |l| we conclude that the QCD suppression of Pq present 
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in the leading order approximation is considerably weakened in the NDR treatment of 
75 after the inclusion of NLO corrections. It is essentially removed for /x > 5 GeV in 
the HV scheme. 

The NLO corrections to Cg which include also the mt-dependent contributions are 
large as seen in table The results in HV and NDR schemes are by more than a 
factor of two larger than the leading order result Cg = Pq^ which consistently should 
not include mt-contributions. This demonstrates very clearly the necessity of NLO 
calculation which allow a consistent inclusion of the important mt-contributions. For 
the same set of parameters the authors of ref. would find Cg to be smaller than 
Cg^^^ by 10-15%. 

The fi and Ajj^ dependences of Cg are quite weak. We also find that the rrit dependence 
of Cg is rather weak. Varying rrit between 150 GeV and 190 GeV changes Cg by at 
most 10%. This weak rrit dependence of Cg originates in the partial cancellation of 



rrit dependences between Y{xt) and Z(xt) in ( 2.10| ) as already seen in the case of 
Kl — > '7i^e~^e~ . Finally, the difference between Cg^^ and Cg^ is small and amounts 
to roughly 5%. 



5.4 p 

5.2 - 

5.0 - 

4.8 - 



„^ 4.4 - 



[1= 2.5 GeV 
\l= 5.0 GeV 
II = 10.0 GeV 



^ 4.2 z 




1.0 
0.9 
0.8 
0.7 



0.6 ^ 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 



0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



Figure 2: Comparison of the Wilson coefficient Cg^^ as a function of s for rrit = 170 GeV, 
^Ms ~ 0.225 GeV and different values of /z in leading order (thin lines) and next-to-leading 



order (thick lines) accuracy. Note the different scales for the real and imaginary parts! 



In fig. I we show Cg^^ of (^^) as a function of s for rrit = 170 GeV, A^g = 225 MeV 



and 2.5 < < 10 GeV. In order to see the importance of the term resulting from the 
one-loop matrix elements one should compare these results with the s-independent values of 
Cg. We should also remember that the NLO corrections to Pq calculated here shift Cg-^'^ for 
fi = 5.0 GeV by ACg^^ ~ 0.8 and ACg^ ~ 0.6 with similar results for other /i. In order 
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to show this effect more exphcitly we also plot in fig. |^ a "leading order" result obtained by 
using only the leading term in (|2.11|) with at the one-loop level but keeping otherwise 
all explicit NLO terms in ( |2.1CI| ) and the contributions from one-loop matrix elements given 
in ( |2.28|) . It should be stressed that roughly 50% of the difference between the "thick" and 
"thin" lines in fig. |^ is due to the term 4/9 in (|3.3|) which in the NDR scheme enters the 



NLO terms in Pq but in the HV scheme is present in the one-loop matrix elements. We have 
left it out in the "thin" lines in fig. ^ in order to show its importance. The calculation of 
NLO corrections to Pq allows a consistent inclusion of this term which contributes positively 
to Cg^^ . Additional enhancement comes from using the two- loop renormalization group 
analysis for Cg and as at the two-loop level. In fig. |^ we also note that ReCg'^'^ ^ ImC'g'^''^. 
The pronounced peak for s = 4m1/m\ = 0.34 is related to the behaviour of h{z, s) in (|2.29|) . 
This peak essentially disappears for n = 2.5 GeV because of the accidential cancellation 
SCi'^ + C2^^ ^ in the dominant term multiplying h{z,s). The authors of ref. [Q] would 
find ReCg^^ by about 15% below our values. In the absence of QCD corrections, h{z,s) in 
]) is multiplied by 6*2°^ = 1 and consequently there is no accidental suppression of this 





term as in the QCD case. Since in addition for as 



tdNDR 
^0 



is slightly enhanced over the 



values given in table |l|, we find Cg^^ in the absence of QCD corrections to be substantially 
larger than the result given in fig. 0. For instance, KeCg^^ varies between 5.2 and 6.3 for 
0.1 < s < 0.9. The complete result for R{s) in this case is shown in fig. ^ at the end of this 
section. 
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a) 



b) 



|X = 2.5 GeV 

11= 5.0 GeV 

II = 10.0 GeV 



- m, = 150 GeV 
- m, = 170 GeV 
m = 190 GeV 



1 

- 12 

- 11 

- 10 

- 9 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 

s 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 

s 



V3> 



Figure 3: a) R{s) for rrit 

b) R{s) for /i = 5 GeV, A- 



170 GeV, A^s 



MS 



225 MeV and different values of 
225 MeV and various values of rrit. 



/i. 



We next present a numerical analysis of ( ^.27| ). In doing this we keep in mind that for 
s ^ m^/m^, s ~ m^i/m^ etc. the spectator model cannot be the full story and additional 
long distance contributions discussed in refs. [|T^|-|T^ have to be taken into account in a 
phenomenological analysis. Similarly we do not include l/^b corrections calculated in 
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which typically enhance the differential rate by about 10%. 

In fig. we show R{s) for rrit = 170 GeV, = 225 MeV and different values of /i. 
In fig. ^) we set = 5 GeV and vary rrit from 150 GeV to 190 GeV. The remaining n 
dependence is rather weak and amounts to at most ±8% in the full range of parameters 
considered. The rrit dependence of R{s) is sizeable. Varying mt between 150 GeV and 190 
GeV changes R{s) by typically 60-65% which in this range of mt corresponds to R{s) ~ m^. 
It is easy to verify that this strong mt dependence originates in the coefficient Cio given in 
( ^.8[ ) as already stressed by several authors in the past 

We do not show the dependence as it is very weak. Typically, changing from 
140 MeV to 310 MeV decreases R{s) by about 5%. 

7 
6 
5 
4 
3 

^ 2 
'o 

Zl 1 

-1 

-2 
-3 
-4 
-5 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

A 

S 

Figure 4: Comparison of the four different contributions to R{s) according to eqn. ( p.27| ). 

R{s) is governed by three coefficients, Cg^^ , Cio and cf''^'^^ . It is of interest to investigate 
the importance of various contributions. To this end we set Aj^ = 225 GeV, mt = 170 GeV 
and /i = 5 GeV. In fig. | we show R{s) keeping only C^" , C^, C^f^^" and the cf^^^^-C^" 
interference term, respectively. Denoting these contributions by -Rg, i?io, Ri and R-j/g we 
observe that the term Rj plays only a minor role in R{s). On the other hand the presence 
of C^"*^^^ cannot be ignored because the interference term Rjig is significant. In fact the 
presence of this large interference term could be used to measure experimentally the relative 
sign of Cj°^'^^ and ReCg^^ which as seen in fig. H is negative in the Standard 

Model. However, the most important contributions are Rg and i?io in the full range of s 
considered. For mt ~ 170 GeV these two contributions are roughly of the same size. Due to 
a strong mt dependence of i?io, this contribution dominates for higher values of mt and is 
less important than Rg for mt < 170 GeV. 
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Figure 5: R{s) for mt = 170 GeV, A^jg = 225 MeV and /i = 5 GeV. 



Next, in fig. || we sfiow R{s) for = 5 GeV, mt = 170 GeV and = 225 MeV compared 
to the case of no QCD corrections and to the resuhs Grinstein et. ah [§] would obtain for 
our set of parameters using their approximate leading order formulae. 

Finally, we would like to address the question of the definition of mt used here. In order to 
be able to analyze this question, one would have to calculate perturbative QCD corrections 
to the functions Y{xt) and Z{xt) and include also an additional order in the renormalization 
group improved perturbative calculation of Pq. The latter would require evaluation of three- 
loop anomalous dimension matrices, which in the near future nobody will attempt. In any 
case, we expect only a small correction to Pq. The uncertainty due to the choice of in 
mt(/i) can be substantial, as stressed in refs. [|33|,0, and may result in 20-30% uncertainties 



in the branching ratios. It can only be reduced if 0{as) corrections to Y{xt) and Z{xt) are 
included. For — > 7r"'"z/i/, Kl — >■ Tr^i/i/, B — and B — >■ X^z/z/ this has been done in 
refs. 



33,34 



The inclusion of these corrections reduces the uncertainty in the corresponding 
branching ratios to a few percent. Fortunately, the result for the corrected function Y{xt) 
given in refs. |33,34| can be directly used here. The message of refs. |]33|,^ is the following: 
For mt = mt(mt), the QCD corrections to Y{xt) and consequently to C*io are below 2%. 
Corresponding corrections to Z{xt) are not known. Fortunately, the mt dependence of Cg is 
much weaker and the uncertainty due to the choice of in mt(/x) is small. On the basis of 

we believe that if mt = mt(mt) is chosen, the 
XoC^e^ ) should be small. 



these arguments and the result of refs. [p^ ,[ 
additional short distance QCD corrections to BR(i? 
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5 Summary 



We have calculated the effective Hamiltonian relevant for the rare decay B Xse^e~ 
beyond the leading logarithmic approximation. The main new result of this paper is the 
calculation of the Wilson coefficient of the operator Qg = (s6)v'_A(ee)y including next- 
to-leading logarithms in the NDR and HV renormalization schemes. A separate analytic 
expression for Cg given in sect. ^ as opposed to C^^^ given in []T2[ should be useful not only 
in 5 — > XgC+e^ but also in i? — > K*e^e~ and other rare -B-decays to which Qg contributes. 
Calculating B X^e+e" in the spectator model we confirm the very recent result for C^^^ 
presented by Misiak in |12|. The cancellation of the scheme dependence in C^^^ is shown 



explicitly in our paper. 

The effect of the NLO corrections is to enhance BR(i? XsC+e^ ) so that its suppression 
found in the leading order analysis of ref. |@| is considerably weakened. This is seen in 
particular in fig. ^. 

We have investigated the mt, and /i ~ 0{m\,) dependence of the "reduced" branching 
ratio R{s). The dependences on A^ and /x are rather small, at most ib8% in the full range 
of parameters considered. The dependence on mt is sizeable. In the range 150 GeV < mt < 
190 GeV it is roughly parametrized by R{s) ~ m^. For mt = 170 GeV, Aj^ = 225 MeV, 
H = 5 GeV and 0.1 < s < 0.8 we find 

1.0 ■ 10"^ < < 9.8 ■ 10"^ (5.14) 

This result can be modified by non-perturbative 1 / corrections and long distance contri- 



butions [|14|-p!q], which are however beyond the scope of this paper. 



Acknowledgment 

We would like to thank Mikolaj Misiak for the correspondence related to his independent 
analysis presented in ref. |T2|. One of us (M.M.) appreciates helpful discussions with M. 
Jamin, M. Lautenbacher and U. Nierste. 



References 

[1] W. S. Hou, R. I. WiLLEY AND A. SONi, Phys. Rev. Lett. 58 (1987) 1608. 

[2] B. Grinstein, M. J. Savage and M. B. Wise, Nucl. Phys. B319 (1989) 271. 

[3] S. Bertolini, F. Borzumati, A. Masiero and G. Ridolfi, Nucl. Phys. B353 (1991) 
591. 

[4] A. Ali, T. Mannel, T. Morozumi, Phys. Lett. B273 (1991) 505. 

[5] W. Jaus and D. Wyler, Phys. Rev. D41 (1990) 3405. 

[6] N. G. Deshpande, K. Panose and J. Trampetic, Phys. Lett. B308 (1993) 322. 

[7] A. Ali, G. F. Giudice and T. Mannel, preprint CERN-TH 7346/94, \hep-ph/ 940821^ . 

[8] C. Greub, a. Ioannissian and D. Wyler, preprint ZU-TH 25/94, \iep-ph/940838i . 



16 



[9] R. GRIGJANIS, p. J. O'DONNELL, M. SUTHERLAND AND H. Navelet, Phys. Lett. B223 
(1989) 239. 

[10] G. Cella, G. Ricciardi and A. Vicere, Phys. Lett B258 (1991) 212. 

[11] M. MISIAK, Nucl. Phys. B393 (1993) 23. 

[12] M. MislAK, Erratum, to appear in Nucl. Phys. 

[13] A. Falk, M. Luke and M. J. Savage, Phys. Rev. D49 (1994) 3367. 

[14] C. S. LiM, T. MOROZUMI and A. I. Sanda, Phys. Lett. B218 (1989) 343. 

[15] N. G. Deshpande, J. Trampetic and K. Panose, Phys. Rev. D39 (1989) 1461. 

[16] P. J. O'DONNELL and H. K. K. Tung, Phys. Rev. D43 (1991) R2067. 

[17] A. J. BuRAS, M. Jamin and M. E. Lautenbacher, Nucl. Phys. B408 (1993) 209. 

[18] A. J. BuRAS, M. MISIAK, M. MUNZ AND S. POKORSKI, Nucl. Phys. B424 (1994) 374. 

[19] A. J. BuRAS, M. E. Lautenbacher, M. Misiak and M. Munz, Nucl. Phys. B423 (1994) 
349. 

[20] G. 'T HooFT AND M. Veltman, Nucl. Phys. B44 (1972) 189. 

[21] M. CiucHiNi, E. Franco, G. Martinelli, L. Reina and L. Silvestrini, Phys. Lett. 
B316 (1993) 127. 

[22] M. CiucHiNi, E. Franco, G. Martinelli and L. Reina, Nucl. Phys. B415 (1994) 403. 

[23] M. CiucHiNi, E. Franco, L. Reina and L. Silvestrini, Nucl. Phys. B421 (1994) 41. 

[24] G. Cella, G. Curci, G. Ricciardi and A. Vicere, Nucl. Phys. B431 (1994) 417. 

[25] G. Cella, G. Curci, G. Ricciardi and A. Vicere, Phys. Lett. B325 (1994) 227. 

[26] N. Cabibbo and L. Maiani, Phys. Lett. B79 (1978) 109. 

[27] C. S. Kim and A. D. Martin, Phys. Lett. B225 (1989) 186. 

[28] M. Jezabek and J. H. Kuhn, Nucl. Phys. B320 (1989) 20. 

[29] R. Fleischer, private communication. 

[30] T. INAMI and C. S. LiM, Prog. Theor. Phys. 65 (1981) 297. 

[31] A. J. BuRAS, M. Jamin, M. E. Lautenbacher and P. H. Weisz, Nucl. Phys. B400 
(1993) 37. 

[32] M. MisiAK, private communication. 

[33] G. BUCHALLA and A. J. BuRAS, Nucl. Phys. B398 (1993) 285. 

[34] G. BUCHALLA and A. J. Buras, Nucl. Phys. B400 (1993) 225. 



17 



